Spontaneous Emission in ultra-cold spin-polarised anisotropic Fermi Seas 
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We examine and explain the spatial emission patterns of ultracold excited fermions in anisotropic 
trapping potentials in the presence of a spin polarised Fermi sea of ground state atoms. Due to 
the Pauli principle, the Fermi sea modifies the available phase space for the recoiling atom and 
thereby modifies its decay rate and the probability of the emitted photon's direction. We show that 
the spatial anisotropies are due to an intricate interplay between Fermi energies and degeneracy 
values of specific energy levels and identify a regime in which the emission will become completely 
directional. Our results are relevant for recent advances in trapping and manipulating cold fermionic 
samples experimentally and give an example of a conceptually new idea for a directional photon 
source. 
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I. INTRODUCTION 

Cold samples of neutral fermionic atoms have become 
an important test-bed for a large number of interesting 
phenomena in many particle physics Since the first 
realisation of quantum degeneracy using two spin com- 
ponents oi'^^K 0, the field has moved quickly from fun- 
damental quantum statistical experiments Q into other 
areas such as BEC-BCS transitionsH , solid state physics 
and even quark-gluon physics [g|7 

A particularly impressive achievement has been the re- 
alisation that by using a Feshbach resonance it is possible 
to incite single atoms to pair 0, [1] . Depending on which 
side of the resonance the experiments are carried out, 
these pairs then either represent bosonic molecules, which 
in turn form a Bose-Einstein condensate, or Cooper- 
pairs, which lead to the formation of a BCS state. By 
sweeping across a Feshbach resonance one can therefore 
explore the BEC-BCS crossover regime, which was until 
recently not experimentally available 

Besides molecular or Cooper pair physics, mono- 
atomic gases have also shown a large potential for demon- 
strating new and exciting physics. While the most dra- 
matic consequence of the antisymmetry condition on the 
wave- function of identical fermions is the formation of the 
Fermi sea at low temperatures, other effects have been 
predicted and been observed. Among them are the mod- 
ification of the scattering properties of two atoms, which 
leads to a reduced efficiency of evaporative cooling [9l.[l0j. 
narrowing of the line-width of light propagating through 
the gas [11 
scattering 
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and the suppression of off-resonant light 
3- 

The inhibition of spontaneous emission in the presence 
of a ground state Fermi sea is another fundamental pre- 
diction which results directly from the Pauli principle 
[U, m, [iBl ■ In it a degenerate Fermi sea of a spin po- 
larised gas forms the environment for a single, excited 
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atom of the same kind. Due to the Pauli principle, the 
Fermi sea effectively blocks out a large amount of the 
phase space that would otherwise be available to the ex- 
cited atom after a de-excitation transition. This leads 
to a modification of the emission properties of the ex- 
cited atom and the details of the effect are determined 
by the size of the Fermi sea, the systems temperature 
and the anisotropy of the trap [3 ■ The influence on the 
lifetime of the excited atom has been recently exhaus- 
tively investigated [3 and the effect was shown to be 
an atom-optical analogue of well known effects in cavity 
QED [13. 

In this work we will investigate the influence of this 
Pauli-blocking effect on the spatial distribution of the 
emission spectrum of a single atom in an anisotropic trap. 
The fact that the emission spectrum becomes anisotropic 
was first shown in and a simple explanation for this 
effect was given. Here we investigate the pattern forma- 
tion in detail and in particular consider highly anisotropic 
traps, which can be achieved today by experimentally us- 
ing, for example, atom chips or optical lattices [l8| . 

In Sec. |TI] we will first describe the model we are using 
and, in Sec. IIIIl derive a relation between the Fermi en- 
ergy and the number of particles in an anisotropic trap. 
In Sec. IIVI we describe the effect of the anisotropy on 
the behaviour of the individual transition elements for 
spontaneous emission and apply the results to explain 
the specifics of the overall emission pattern. Finally we 
conclude. 



II. MODEL 

We consider an ideal gas of spin polarised fermions 
trapped in a harmonic potential. All atoms are assumed 
to be in their internal ground state, \g), so that the gas 
becomes quantum degenerate at low enough tempera- 
tures and forms a perfect Fermi sea at absolute zero. 
In the following we will restrict our calculations to this 
limit, as in it the effects we describe are most pronounced 
and the extension to finite temperatures is, while com- 
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putationally challenging, conceptually straightforward. 

In addition to the Fermi sea we assume the presence 
of a single extra fermion, which is distinguished from the 
others by being in an internally excited state, |e). After 
some time this atom will spontaneously emit a photon, 
make a transition into the ground state and become part 
of the Fermi sea. As all atoms are assumed to be spin po- 
larised, the Pauli principle demands that the new ground 
state atom has to join the Fermi sea with an energy larger 
than the Fermi energy. This is an energetically very un- 
favourable process and the presence of the Fermi sea will 
therefore lead to an inhibition of the spontaneous emis- 
sion rate with respect to the case of a free space particle 

In the following we will denote the spontaneous emis- 
sion rate of photons along the direction O and into the 
solid angle dfl in the presence of ground-state fermions 
by r{n) dn and compare it to the free case {N = 0) which 
we denote by To{fl)d^. Using Fermi's golden rule we can 
express the excited atom's decay rate as 



J2 P,n{l ~ F^)\{n\e-^^'^''y^^rn)\\ (1) 



n,m— 



where F„ = (^eif^^/kBT){\-n) ^ ^^-i ^^le Fermi-Dirac 
distribution function and (1 — i^„) is the probability that 
an energy level \n) of the harmonic trap is unoccupied. 
P„ = p^f,~itiuj/kBT)X-rh jg ^Yie Bohzmann distribution 
function describing the single excited fermion in state 
m) of the harmonic trap, which in turn is assumed to 
have the frequencies (wa; , a;^ , a;^ ) = ujX. If we restrict 
eq. ([1]) to zero temperature the Fermi-Dirac distribution 
function becomes a step function and hence only states 
with an energy greater than the Fermi energy have a fi- 
nite value for 1 — F^- Similarly, the excited fermion will 
occupy the ground state of the harmonic trap, \m) — |0), 
and eq. ((T|) simplifies to 
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where np represents the Fermi shell. The inhibition of 
spontaneous emission that results from equation ([T]) has 
recently been investigated [l6| and here we extend this 
work by presenting a thorough and detailed investiga- 
tion into the spatial emission patterns that result from 
it. This is of interest due to the new parameter ranges 
which have become experimentally available in recent 
years. These include, in particular, highly anisotropic 
traps. While the appearance of a fine structure in the 
emission patterns was already mentioned in , here we 
derive the framework for its explanation. 

We assume that the harmonic oscillator potential has 
the following standard form 



V 



{X-rf 



where M is the mass of the particle and the values of 
Xx,y.z determine the degree of anisotropy in the different 
directions. For numerical simplicity we will only deal 
with values of A > 1 and restrict ourselves to two types 
of anisotropic trapping potentials in which two of the 
axes have identical strength. 



Xx = Xy = I and Xz = X, 



X-, 



Xy = X and Xz 



pancake shape , (4a) 
cigar shape . (4b) 



Due to the symmetries of the pancake and cigar shaped 
traps about the tight and the soft axes, respectively, 
eq. (m can be simplified and expressed in terms of in- 
complete gamma functions for both trap shapes 



Mfi9) 



7(nF + l,/3) 
r(nF + l) 



n=0 



/9" 7(L^J+l,f) 

r(L^J + i) 



(5) 

Here [xj denotes the largest integer less than or equal 
to X and ni;- is the quantum number of the Fermi shell, 
which we will describe in more detail later. For brevity 
we have defined 

a ~ if' cos^ 6*, P — svp? 9 pancake shape , (6a) 
a = if sin^ 6*, P ^ nf cos^ 9 cigar shape , (6b) 

in which r/ represents the Lamb-Dicke parameter. One 
can immediately see that the angular distributions for 
the pancake and the cigar shaped trap can be obtained 
from each other by a simple 7r/2 rotation, as one would 
expect. We will make use of this fact when discussing 
emission patterns in Sec. IIVI 



III. DEGENERACIES 

Let us first discuss the relationships between the 
different parameters characterising a Fermi-sea in an 
anisotropic trap. Since the degeneracies of states with 
equal energy are a function of the trapping frequencies 
in the different directions, the relationship between the 
Fermi energy and particle number is not as straightfor- 
ward as in the well-known isotropic case. 

The eigenenergies of the harmonic potential in eq. ([3]) 
are given by 
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(7a) 
(7b) 



[Xix' + Xly' + Xiz'), (3) 



where we have defined the shell quantum numbers of the 
pancake and cigar shaped harmonic traps as n-p — Ux + 
Uy + Xuz and ric = Xux + Xny + Uz, respectively. As usual, 
nx,ny and n, refer to the integer quantum numbers of 
the harmonic oscillator. 

As the aspect ratio of a trapping potential is increased 
the resulting energy levels typically have a reduced de- 
generacy relative to the isotropic case, A = 1 [l^. For the 
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purposes of this work, and without loss of generahty, we 
will consider only integer values of A, allowing us in turn 
to restrict ourselves to integer values for rip.c- Therefore 
we can write the degeneracy for states with fixed energy 
as 

ffrip = ^(^p + l)(2np - Arip + 2) , (8a) 



5nc = 2(^c + l)(?^c + 2) , 



(8b) 



where here and in the following all quantities carrying 
a tilde take the value of the quantity without the tilde 
divided by A and rounded down to the nearest integer, i.e. 
X = [a;/AJ. Consequently, the total number of quantum 
states with an energy equal to and smaller than Enp is 
then given by the sum. 



which can be calculated to be given by 
Sp = ^{np + l){2nF - Xnp + "2) 



(9) 
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X^np[2 + hp) 
8 + Snp - AXhp 



+ 3 (10a) 



Sc = -{fip + l){np + 2){'Snp -2fipX + 'S). (10b) 

In our model we assume a spin polarised gas in which each 
oscillator state is filled with one fermion only. Eqs. ([9|) 
therefore determine the number of particles confined for 
a given Fermi energy Ep ~ nphuj + Eq, where Eq is the 
ground state energy of the potential. 



IV. EMISSION PATTERNS IN ANISOTROPIC 
TRAPS 

A. Emission Probabilities. 

To understand the emission patterns later on, let us 
first have a brief look at the emission probabilities of 
an excited atom in an anisotropic trap. In the pres- 
ence of an anisotropic Fermi sea the rate of spontaneous 
emission along a specific direction is determined by three 
parameters: (1) the number of ground state atoms, (2) 
the degeneracy of the available states and (3) the Lamb- 
Dicke parameter 77 = Eji/hcu. The latter determines 
the range of accessible states and is given by the ratio 
between the recoil energy, Ej^ = h^k'^/2M , and the trap- 
ping strength, hjjJx,y,z, in the different directions. Here 
fco is the wave vector corresponding to the transition 
|e) - \9)- 

In this section we will focus on the influence of the de- 
generacies of the available states and therefore on the 
anisotropy of the trap. Let us do this by examining 
the matrix elements for individual transitions from the 
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FIG. 1: Top row: emission probability, Pe, into individual 
shells in a pancake shaped trap for rf — 36 and A = 10, 23, 46. 
The arrow indicates the n = 20 energy level of the harmonic 
trap, which is referred to in the text. Bottom row: decay rate 
of the excited particle, Mf, for the same parameters as above. 



0.06 
P 



0.03 



0.5 



X= 10 








■ X^A6 












25 50 75 „ 


25 50 75 „ 


25 50 75 „ 


\ 




V 







25 50 75 „ 25 50 75 „ 25 50 75 „ 

F F F 



FIG. 2: Top row: emission probability, Pe, into individual 
shells in a cigar shaped trap for rj^ — 36 and A = 10, 23, 46. 
The arrow indicates the n = 20 energy level of the harmonic 
trap, which is referred to in the text. Bottom row: decay rate 
of the excited particle, Mf, for the same parameters as above. 



ground cm-state of the excited atom to a single final 
state, |n). 



ikx I n\ |2 



(11) 



It is well known that for an isotropic trap this distribu- 
tion is Gaussian in shape and centered around an energy 
level n = rf' . The effects introduced by an anisotropy are 
significant and can be clearly seen in the graphs in the 
upper rows of Figs. [T] and [2l where we show for a pan- 
cake and a cigar shaped trap, respectively, for increasing 
values of the anisotropy, A = 10, 23 and 46. The most 
obvious feature in both situations is the appearance of 
a A-dcpcndent discontinuity in the distribution, which is 
more pronounced in the cigar shaped setting. 

To explain this behaviour, let us first intuitively argue 
its existence. When an internally excited atom which is 
trapped in the ground state of an empty isotropic har- 
monic trap decays, the probability of the photon being 
emitted is the same in all directions. This is rather easy 
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FIG. 3: (a) Emission probability into individual states within 
the shell n = 5 in an isotropic trap. The triplets represent 
{nx,ny, riz) and all permutations of each triplet have the same 
probability. The Lamb-Dicke factor is 77 = 5, however chang- 
ing this value only scales all values, (b) Emission probability 
into individual states within the shell n = 5 in a pancake trap 
with A = 5. The values for the three states on the left are not 
visible on this scale. 



probability of occupation (see Fig. [3Kb)) the overall emis- 
sion probability into this energy shell is increased, leading 
to the observed discontinuous jump. As an example let 
us consider a pancake shaped trap with an aspect ratio of 
A = 10. For the shell n = 19 (indicated by the arrow in 
Figs. [T]) the degenerate states are {rix + Uy, n^) — (19, 0) 
and (9,1), whereas for the n = 20 energy level states 
are (n^ + ny,n^) = (20, 0), (10, 1) and (0,2). The ex- 
tra (0, 2) state is the dominant contributor and its ap- 
pearance responsible for the discontinuous increase in the 
emission probability. For the cigar trap this effect is even 
more pronounced as there are two tight directions and in 
the example above the states {nx,ny,nz) = (2,0,0) and 
(0,2,0) become both available. 

For completeness we show the integrated emission 
probability for increasing particle number (i.e. increasing 
Fermi energy or Fermi level) and different anisotropies in 
bottom rows of Figs. [1] and [2l Fermi inhibition is absent 
for the empty trap (Mf = 1), shown for np = —1, slowly 
increases for n^? > and accelerates for np ^ rj^ . The 
discontinuity in the variable Pe {n) translates clearly into 
non-smooth kinks in this distribution. 



B. Emission along Tight and Soft Axes. 



to understand as in this situation the density of states is 
identical in all directions. However, for the anisotropic 
trap the situation is different. As the aspect ratio is in- 
creased, the degeneracy of any energy level will either 
decrease or remain the same. Therefore, up to a specific 
shell the number of quantum states, given by eqs. ([9]), 
will be reduced and, as a result, the density of states in 
the different directions changes. Since the recoil of the 
de-excited fermion to a certain quantum state and the 
direction of the emitted photon are directly related, it 
seems rather surprising that for the free case this emission 
remains isotropic irrespective of the diminishing number 
of quantum states. However, it is exactly the modified 
distribution shown in the upper rows in Figs. [T] and [2] 
that makes this phenomenon possible. 

To gain more insight into the source of the discontinu- 
ities let us consider the emission probability into specific 
states within a degenerate shell n of an isotropic and a 
pancake shaped (A = 5) trap. Fig.[3]shows Pe{nx, riy, Uz) 
for a fixed shell, in which all combinations of the triplet 
(in both (a) and (b)) of quantum numbers adds up to 
n = 5. It can be seen that in general states which in- 
clude ground state excitations have a higher probability 
for occupation than the ones which do not, which is due 
to the fact that the excited atom is initially in its centre- 
of-mass ground state. 

When we move from the isotropic to the anisotropic 
setting it is therefore clear that whenever the value of 
an energy shell, n, reaches an integer multiple of the 
anisotropy parameter, the shell contains a state with two 
ground state excitations. As these states have a higher 



The fact that the presence of an anisotropic Fermi 
sea will lead to anisotropic emission patters was already 
noted in and in the following we will develop a de- 
tailed understanding of the directional features. Since 
for the pancake as well as for the cigar shaped trap the 
emission is isotropic around their respective symmetry 
axis, (0, tt), we can treat both geometries in a quasi 2D 
picture. It is then immediately clear that the results for 
both settings will be related by a simple 7r/2 rotation 
(due to our definition of A > 1). 

Let us first look at the emission along the principal 
axes of the anisotropic trap in the tight and the soft 
direction. Choosing the tight direction in the pancake 
(cigar) shaped trap along 6 — {6 ~ tt/2) the modifica- 
tion factor in eq. ^ simplifies to 
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(12) 



The behaviour of this equation with increasing 
anisotropy is shown in Fig. |4] for a system with np = 
60. The most obvious feature of the plot is a scries of 
sawtooth- like discontinuities. Careful examination shows 
that n^^ of these exist and they appear whenever the value 
of the aspect ratio. A, increases beyond the values of — , 
(jn = 1,2, ...,ni?). The increase in emission probability 
for values just after this point is due to the availability 
of an extra free state with a lower tight excitation just 
outside the Fermi edge. For example, in the pancake 
trap, when one moves from A = 30 to A > 30 the state 
{rix -\- nyjfiz) = (0,2) emerges from the Fermi sea for 
np — 60. As discussed in Section flV Al this state has a 
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FIG. 5: (a) Mf{6) in a pancake shaped trap at T = 0. = 
25, A = 4 with np ~ 31 (outermost) to np ~ 36 (innermost). 



(b) Mf{9) in a pancake shaped trap with A 
and np = 23. 
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FIG. 4: M/ along the tight axis at T = 0. rf = 49, np = 60. 
Note that we use a continuous distribution of A for this graph. 



high probability to be emitted into as it contains ground 
state excitations in the soft direction, hence the large in- 
crease in the decay rate. By increasing A further this 
state moves away from the Fermi edge, and the emission 
probability decreases until the next state with a lower 
tight excitation, emerges from the Fermi sea. For val- 
ues of A > np no more discontinuities appear since the 
Fermi sea only occupies energy states with ground state 
excitations in the tight direction. Emission along the 
soft direction can be calculated from eq. ([12]) by taking 
A = 1 and the decay rate along this direction is deter- 
mined exclusively by the Fermi shell np and the value of 
the Lamb-Dicke parameter rj. 

Considering a fixed value of the aspect ratio A in ei- 
ther anisotropic trapping potential and changing np one 
notices a degeneracy in the emission probability in the 
tight direction, shown in Fig. [5]Ja). This behaviour was 
already mentioned in and we can see from eq. (jl2p 
that it stems from the fact that fi p only changes its value 
in steps of A. An increase in the value of coincides 
with the Fermi sea occupying a state with a higher tight 
excitation (and ground state soft excitations), leading to 
a decrease in the decay rate along the tight direction. 
For example, when moving from np = 35 to rii? = 36 
the state (n^ -I- n^, n^) = (0, 9) becomes occupied by the 
Fermi sea, producing the discontinuous reduction of the 
decay rate (see Fig. [S^a)). 



C. Fine Structure 

The emission spectrum between the principal axes is 
characterised by the appearance of a fine structure (see 
Fig. Efb)), which exists for a wide range of parameters. 
The first hint to understanding the origins of the visible 
extrema comes from noticing that the number of maxima 
between the soft and tight axes is related to the number 
of excitations in the tight direction that are occupied by 



the Fermi sea, fip. To show this relation let us consider 
the emission probability into shells with a fixed value for 
nz in a pancake shaped trap 
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nzl r(max(0, np — Xnz + 1)) 

(13) 

with the definition of a and /3 given in eqs. ([6]). As a spe- 
cific example we show in Fig. \5lh) a gas with np = 23, 
rj^ = 25 and A = 11. In this case we find np = 2 maxima 
in the tt/2 arc between the tight to the soft axis. Com- 
paring this emission pattern to the results from eq. p3p . 
one can see (Fig. [6|) that each isolated contribution from 
a transition into a state with a fixed value of is respon- 
sible for one of the maxima. For values of > hp the 
emission is predominantly into the tight direction, there- 
fore originating from transitions into states for which 
both ground state excitations in the soft direction are 
available. Similarly, when restricting the recoiling atom 
to occupying states with a ground state excitation in the 
tight direction, nz = 0, the emission is mainly focussed 
around small angles about the soft axis. The intermedi- 
ate excitations, = 1, 2, make up the two intermediate 
ripples between the principle axes and summing up the 
contributions to the photon emission of all four plots in 
Fig. [5] gives the emission plot shown in Fig.[5]^b). In con- 
trast, if we calculate eq. for an isotropic trap for 
different values of , each individual term would show a 
similar behaviour of having a single maximum at a finite 
angle between the principle axes. However, the sum of 
those will give the isotropic emission pattern which corre- 
sponds to the decay rate being the same in all directions. 

It is now obvious that for the limit \ > np the fine 
structure disappears and the extrema of emission will be 
located around the directions of the principal axes (see 
Fig. EJa)). As A — > 00, emission into the tight direction 
is reduced, whereas the emission in the soft direction re- 
mains constant. Fig. [7]^b). In this regime the Fermi sea 
is completely confined to states with ground state excita- 
tions in the tight direction. Therefore, it becomes easier 
for the recoiling atom to access states in the soft direc- 
tion due to the diminishing density of states in the tight 
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FIG. 7: (a) M/(0) in a cigar shaped trap at T = with 
np = 45 and r]^ = 49. A = 46 (outermost), A = 96 (center 
plot) and A = oo (innermost). The plot is symmetric through 
a 2tt rotation about the (0, tt) axis, (b) A three-dimensional 
illustration of the excited particles decay rate in a cigar trap 
in the large anisotropy limit. 
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FIG. 6: Mf{0,nz) for a pancake shaped trap with, A = 11 
rj^ — 25 and np = 23. In the four graphs the decay is only 
allowed into quantum states of the harmonic trap with riz = 
0, 1, 2, and Uz > 3, respectively. 



direction. In the limit of A 
can be written as 



oo the emission probability 



M/(6l; A -> oo) = 



r(nF + l) 



(14) 



and it shows that the emission probability in the tight 
direction has completely vanished. 

It is possible to make use of this behaviour and cre- 
ate a system where photon emission become highly di- 
rectional. While directional photon emission is usually 
achieved by using optical cavities (and therefore engi- 
neering the Hilbert space of the photon), this example 
is complementary in that it uses a cavity (trap) for the 
atoms and thereby engineers the Hilbert space of the par- 
ticles. Let us stress that it is not primarily the size of the 
Fermi sea that is responsible for this effect, mearly the 
presence of the Fermi sea. The emission probability of 
the photon can still be close to the emission probability 
in free space whilst 77^ > np. (see Fig. [Tl^a)). As the 
emission is symmetric through a 2tt rotation about the 
(0, tt) axis in the above example, we display the 3D emis- 
sion probability in Fig.[7jb). Also note that for a pancake 
shaped trap this effect would correspond to emission into 
a well defined plane perpendicular to the tight principal 
axis. 



CONCLUSION 



In this work we have given a detailed investigation into 
the spatial properties of spontaneous emission of a single 
atom in the presence of an anisotropic, ideal and spin 
polarised ultracold Fermi gas. The demand to obey the 
Pauli principle leads to the formation of a non-trivial, 
anisotropic emission pattern for the photon, which can be 
explained by carefully examining the allowed transitions 
the recoiling atom can make. 

We have first calculated the relation between the Fermi 
energy and the particle number and then investigated the 
single particle transition matrix element, for both geome- 
tries of anisotropic traps. The change in the density of 
states into the different spatial directions was found to be 
accompanied by the appearance of discontinuities in the 
distribution of the emission probability spectrum for dif- 
ferent shells. While in an isotropic trap these two effects 
cancel and produce an isotropic emission spectrum, in an 
anisotropic trap they lead to an intricate fine-structure 
in the presence of a Fermi sea. 

In a next step we have managed to explain this fine- 
structure by attributing the extrema to the emissions 
which come from the transitions of the recoiling atom into 
well defined states in the tight direction. If the aspect ra- 
tio exceeds the Fermi energy, the fine-structure vanishes 
and the emission spectrum becomes smooth, though not 
isotropic, again. 

Finally, we have pointed out that this system can be 
used to create a highly directional photon source. The 
effect uncovered is complementary to the common use 
of optical cavities to influence a photons direction after 
emission and makes use of the ability to influence the 
atom's phase space. The experimental observation of di- 
rectional photon emission in anisotropic, cold, fermionic 
gases would therefore be a sign of a fundamental conse- 
quence of the symmetry of fermionic particles. 
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